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ABSTRACT 


The  effects  of  multipathing  on  the  Rayleigh  spectrum 
mask  depth-dependent  spectral  zeroes  and  cause  significant 


errors  in  magnitude  determinations.  A  new  technique,  called  ..  jj 
phase  equilization  filtering,  is  introduced  for  the  purpose  -I 

of  eliminating  multipath  effects.  The  technique  is  demonstrated 
to  be  effective  on  synthetic  signals  in  which  the  undisturbed 

spectrum  is  known.  In  a  typical  ca so  an  18-dB  spectral  hole 

. 

at  a  period  of  20  seconds,  caused  by  multipathed  energy  arriving 
only  30  seconds  later  than  the  primary  arrival,  was  essentially 
eliminated.  The  usefulness  of  the  technique  is,  also  demonstrated 


by  application  to  teleseismic  Rayleigh  waves.  . "r  . ' 
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Introduction 


The  presence  of  non-horizontal  boundaries  in  the  crust 
and  upper  mantle  of  the  earth  result  in  reflection  and  refrac¬ 
tion  of  Rayleigh  and  Love  waves.  Secondary  arrivals,  or  multi¬ 
paths,  are  produced,  and  these  cause  constructive  interference 
at  some  frequencies  in  the  recorded  signal  and  destructive 
interference  at  other  frequencies.  The  effect  of  multipathing 
on  the  surface  wave  is  sufficiently  great  to  mask  or  distort 
signal  characteristics  which  are  important  in  the  analysis  of 
seismic  events,  especially  those  originating  in  the  structurally 
complicated  Sino-Soviet  sub-continent.  Multipathing  and  related 
effects  of  the  propagation  path,  unless  identified  and  removed, 
will  limit  the  success  of  research  efforts  directed  toward  iden¬ 
tifying  source  differences  between  explosions  and  earthquakes, 
in  utilizing  the  improved  horizontal  data  expected  from  borehole 
long  period  instruments,  and  in  achieving  accurate  yield  deter¬ 
minations  from  recordings  of  surface  waves. 

In  this  report  we  discuss  a  new  approach  to  the  removal 
of  the  effects  of  multipathing  from  surface  wave  signals.  We 
propose  to  develop  and  to  apply  this  technique,  called  phase 
equilizstion  filtering,  to  Rayleigh  and  Love  waves  originating 


from  events  in  the  Asian  continent. 
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Phase  Equalization  Filtering 


If  the  phase  of  a  dispersed  wave  train  is  known  or  can 
be  determined,  a  filter  can  be  developed  which  has  the  same 
phase  spectrum  as  the  dispersed  signal  and  any  desired  ampli¬ 
tude  spectrum.  Cross-correlation  of  the  signal  and  the  phase 
equalized  filter  produces  a  function,  called  a  pseudo-autocorre¬ 
lation  function  (PAF),  which  possesses  several  desirable  pro¬ 
perties.  The  PAF  contains  the  spectrum  of  the  signal  in  a 
retrievable  form.  Time  compression  of  the  signal  produces  a 
S/N  improvement  on  the  order  of  10  dB  and  allows  better  separ¬ 
ation  of  the  primary  and  multipaths  or  second  signals.  Mathe¬ 
matically,  the  PAF  is  an  even  function  and  has  many  of  the 
desirable  properties  of  the  auto-correlation  function. 

Let  c (i)  be  the  filter  obtained  from  an  estimated  disper¬ 
sion  curve  and  3«  the  recorded  signal.  Then 

cflftgc#  =  £g(t)'c(t+T)dt  =  H(T) 

where  m  is  the  cross-correlation  function.  The  Fourier  trans¬ 
form  of  H(t)  can  be  written 

Q  (w)  •  C  (w)  sxp  i  [&j  M  -  0C  (*>)] 
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Clearly  the  real  part  of  the  transform  of  H(r)  will  be  maximized 
when  the  phase  spectra  of  the  filter  and  the  signal  are  equal. 

In  this  case  the  spectrum  &C")C(w)co3(o)  will  be  real  and  even, 
and  H(j) is  defined  to  be  a  PAF.  The  process  of  maximizing 
the  integral  of  the  real  spectrum  corresponds  in  the  time  domain 
to  maximizing  the  peak  value  of  the  correlation  function.  If 
no  multipathing  is  present,  this  maximization  can  be  achieved 
in  practice  by  varying  the  phase  of  the  filter  unti 1  it  matches 
that  of  the  signal.  If  get)  consists  of  a  primary  and  secondary 
signals  arriving  at  different  times,  HCt)  will  consist  of  a  series 
of  peaks  corresponding  to  each  component  of  the  total  signal, 
the  first  peak  corresponding  to  the  primary  arrival.  Theore¬ 
tically,  Hfr)  endures  from-® to  *00.  In  practice,  for  band-limited 
data  (periods  of  15-70  seconds,  for  example),  the  useful  in¬ 
formation  is  almost  entirely  contained  in  a  window  +  150  seconds 
from  the  peak  of  the  correlation  function.  A  300-second  separ¬ 
ation  in  the  arrival  times  of  a  primary  and  a  mult. i  path  or  be¬ 
tween  different  signals  would  result  in  noninter fern ng  peaks 
and  would  allow  the  determination  of  the  phase  equalized  filter 
by  maximizing  the  amplitude  of  the  primary  peak.  Simultaneously, 
the  spectrum  of  the  primary  signal,  undistorted  by  multipathing, 
can  be  obtained  by  transforming  a  300-second  window  centered 
at  the  primary  peak.  Indeed,  this  process  can  be  accomplished 
not  only  for  the  primary  arrival,  but  for  any  multipath  or 
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interfering  secondary  component  making  up  the  total  signal. 

As  time  separation  between  the  primary  and  secondary 
arrivals  decreases  to  less  than  300  seconds,  use  can  be  made 
of  the  property  of  PAFs  that  in  the  absence  of  multipathing 
they  are  even  functions.  Thus,  in  the  presence  of  multipathing, 
we  know  that  the  right-hand  side  of  the  primary  peak  which  is 
disturbed  by  multipathing  would,  if  the  multipathing  were  not 
present,  be  the  mirror  image  of  the  left-hand  side  which  is 
not  disturbed  by  multipathing.  Replacement  of  the  right-hand 
side  by  the  reflected  left-hand  side  of  H(r)  removes  the  effect 
of  the  multipath,  unless  the  time  separation  between  the  arrivals 
is  so  small  that  the  multipath  peak  affects  the  left-hand  side 
of  the  primary  peak.  In  practice,  a  time  separation  of  about 
50  seconds  appears  to  be  resolvable  for  signals  containing 
periods  up  to  70  seconds. 

An  earthquake  occuring  in  Tadzhik,  USSR,  and  recorded 
at  the  Large  Aperture  Seismic  Array  ( LAS A )  in  Montana  is  use¬ 
ful  in  demonstrating  both  the  profound  effects  of  multipathing 
and  the  effectiveness  of  phase  equalization  filtering  in  re¬ 
moving  the  multipaths.  Figure  1  shows  the  vertical  long-period 
seismograms  recorded  at  sites  AO,  FI  and  F3,  respectively,  at 
LASA,  The  average  properties  of  the  propagation  paths  from  an 
earthquake  10,000  kilometers  away  to  each  of  these  sites,  which 
are  separated  by  about  100  kilometers,  cannot  be  very  different? 
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however,  the  seismograms  are  obviously  different.  The  amplitude 
spectra  of  the  seismograms  recorded  at  AO,  FI,  and  F3,  shown 
in  Figure  2,  are  also  different.  Using  the  criteria  developed 
above,  it  was  possible  to  find  a  single  phase  equalized  filter 
which  fit  each  seismogram.  Correlation  of  this  filter  with 
each  of  the  three  seismograms,  replacing  the  right-hand  sides 
of  the  primary  peaks  of  the  resulting  correlation  functions  with 
the  reflected  left-hand  sides  to  remove  the  multipaths,  taking 
the  Fourier  transform  of  the  PAFs,  and  dividing  the  resulting 
spectra  by  the  spectrum  of  the  filter  resulted  in  the  spectra 
shown  in  Figure  3.  Note  that  these  spectra  with  the  multipath 
effects  removed  are  now  very  similar,  as  might  be  expected  for 
seismograms  representing  almost  the  same  propagation  path.  A 
very  slight  difference  in  propagation  path  can,  and  often  does, 
result  in  dramatic  differences  in  multipath  effects. 

Another  example  from  LASA  demonstrates  the  ??/N  and  time 
compression  aspects  of  phase  equilization  filtering.  A  small 
magnitude  K-i'O  earthquake  whose  epicenter  was  located  approx¬ 
imately  1400  kilometers  northeast  of  the  Tadzhik  event  is  barely 
recognizable  on  the  AO  vertical  LP  channel,  shown  in  Figure  4. 
Figure  4  also  shows  the  function  obtained  by  correlating  the 
signal  with  the  same  filter  determined  for  the  Tadzhik  event. 

The  increase  in  S/N  is  apparent.  Both  the  least  time  arrival 
and  a  multipath  or  aftershock  arriving  140  seconds  later  are 
now  clearly  visible. 
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Phase  equalization  filtering  can  be  used  to  separate 
multiple  events  as  well  as  primary  and  multipath  components. 

A  simulated  mixed  signal  was  obtained  by  adding  the  Tadzhik 
signal  as  recorded  at  site  AO  at  LASA  to  itself  delayed  by 
700  seconds  and  reduced  in  amplitude  by  12  dB .  Figure  5 
shows  the  real  signal,  the  delayed  and  diminished  version  of 
the  same  signal,  and  the  sum  of  the  two.  Note  that  the  presence 
of  the  delayed  event  in  the  coda  of  the  signals  is  not  at  all 
obvious.  Using  the  same  filter  that  was  previously  obtained 
for  the  Tadzhik  signal,  correlation  functions  were  computed  d 

for  both  the  actual  and  the  combined  signal.  These  are  shown  J 

42 

■j£ 

in  Figure  6.  The  primary  signal  and  a  multipath  with  a  delay  | 
of  500  seconds  are  obvious  in  the  real  signal  correlation  function; 
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and  the  peak  corresponding  to  the  simulated  signal  at  a  delay 
of  700  seconds  is  also  obvious  in  the  combined  signal  correla¬ 
tion  function.  Time  separation  of  the  components  of  the  com¬ 
bined  signal  is  so  well  accomplished  that  it  is  now  possible,  ~ 

by  centering  300-secona  windows  on  the  respective  correlation 
peaks,  to  obtain  the  amplitude  spectrum  of  the  primary  signal  ' 

undistorted  by  either  the  multipath  or  the  second  arrival,  and  | 
also  the  amplitude  spectrum  of  the  simulated  second  arrival  I 

undistorted  by  either  the  multipath  or  the  main  signal  coda.  J 

These  two  spectra  are  shown  in  Figure  7.  The  spectrum  of  the  j 
second  arrival  is  properly  12  dB  less  than  that  of  the  main  signal! 
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AMPLITUDE  SPECTRUM  OF  AO  TADZHIK 
WITH  200  SEC  WINDOW  CENTERED 

AT  MAIN  PEAK - 

2nd  PEAK  - 
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Even  the  spectrum  of  the  multipath,  shown  in  Figure  3,  is 
easily  obtainable  because  of  the  time  compression.  The  iso¬ 
lation  and  analysis  of  multipath  spectra  may  make  possible 
'horizontal  exploration  seismology',  the  study  of  the  location 
and  properties  of  lateral  boundaries  in  the  crust:  and  upper 
mantle . 

It  is  constructive  to  consider  the  use  of  phase  equiliza- 
tion  filtering  on  synthetic  signals.  In  this  case  the  signal 
spectrum  prior  to  multipath  distortion  is  exactly  known  and 
can  be  directly  compared  to  the  filtered  estimate.  Consider 
the  synthetic  signal  shown  in  the  top  of  Figure  9.  If  the 
amplitude  of  the  first  600  seconds  of  this  signal  is  zeroed 
and  the  remainder  is  delayed  30  seconds  and  added  to  the  orig¬ 
inal  the  'multipathed '  signal  shown  in  the  bottom  of  Figure  9 
is  produced.  This  record  simulates  a  situation  .in  which  only 
signal  periods  less  than  about  30  seconds  have  been  multipathed 
and  these  have  travelled  less  than  l,k  farther  than  the  primary 
arrivals.  The  effect  of  the  multipath  is  seen  to  double  the 
amplitudes  of  periods  around  30  seconds  and  almost  to  cancel 
periods  around  20  seconds.  Figure  lo  shows  the  ainptitude 
spectra  of  the  two  seismograms.  A  6  dn  'build-up'  at  30  second 
a  10  dB  'build-up'  at  15  seconds,  and  a  1.6  dD  'hole'  at  20 
seconds  are  obvious.  A  magnitude  determination  on  this  signal 
would  be  misleading  because  of  the  mult ipathiny . 
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Because  of  the  small  delay  (30  seconds)  relative  to  the 
primary  with  which  the  multipath  arrives,  it  is  beyond  the 
resolution  of  standard  techniques  to  separate  the  two.  For 
example,  a  multiple  filter  analysis  of  the  simulated  multi- 
pathed  signal  is  shown  in  Figure  11.  Not  only  are  the  multi- 
path-induced  build-ups  and  hole  present,  there  is  no  hint  that 
they  are  caused  by  a  multipath. 

Figure  12  shows  the  function  obtained  by  correlation  of 
the  multipathed  signal  and  the  phase  equalized  filter;  the  filter 


for  the  simulated  signal  of  course  is  known  exactly.  Also 
shown  is  the  PAF.  i.e.,  the  correlation  function  with  the  right 


hand  side  replaced  by  the  reflected  left-hand  side.  The  spectrum 
of  the  PAF  normalized  for  the  filter  spectrum  is  shown  in  Figure 
13,  with  the  spectrum  of  the  multipathed  signal  also  presented 
for  comparison.  The  build-ups  at  15  and  30  seconds  and  the  hole 


at  20  seconds  have  been  strongly  attenuated.  |j 

To  demonstrate  that  phase  equalization  filtering  will  rj 

retain  spectral  holes  which  might  be  present  prior  to  multipathing,  | 


it  is  useful  to  consider  another  simulated  situation.  Figure  14 
shows  a  synthetic  signal  with  a  time-domain  modulation  which 
will  produce  a  spectral  hole.  The  hole  can  be  thought  of  as 
being  produced  by  the  source.  If  the  amplitude  of  the  first 
100  seconds  of  this  signal  is  zeroed  and  the  remainder  is  de¬ 
layed  70  seconds  and  added  to  the  original,  the  'multipathed' 


signal  shown  in  the  bottom  of  Figure  14  is  produced.  The 
spectra  of  the  two  seismograms  are  shown  in  Figure  15.  The 
source  hole  at  a  period  of  22  seconds  is  relatively  unaffected 
by  the  multipath,  but  another  significant  hole  is  introduced 
at  a  period  of  about  50  seconds.  A  correlation  function  is 
obtained  by  cross-correlating  the  mu 1 tips the d  signal  and  the 
phase  equalized  filter.  This  function  and  the  corresponding 
PAF  are  presented  in  Figure  16.  The  spectrum  of  the  PAF,  along 
with  the  spectrum  of  the  multipat.hed  signal,  is  shown  in  Figure 
17.  Nearly  all  of  the  multipath  hole  has  been  eliminated  from 
the  PAF  spectrum;  nearly  all  of  the  source  hole  has  been  re- 
tained,  and  the  period  at  which  the  source  hole  occurs  has 
been  distorted  by  less  than  a  second. 

These  preliminary  experiments  using  both  real  and  synthetic 
data  indicate  that  phase  equalization  filtering  has  great  po¬ 
tential  in  removing  effects  which  are  undoubtedly  causing  scatter 
and  error  in  Rayleigh  wave  discrimination  criteria  and  which  are 
masking  spectral  manifestations  of  the  source.  The  technique 
is  completely  applicable  to  Love  waves  and  can  bo  applied  as 
high-quality,  horizontal  data  becomes  available. 


